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An analytic expression for the '^5'o phase shifts in nucleon-nucleon scattering is derived in 
the context of the Schrodinger equation in configuration space with a short distance cutoff 
and with a consistent power counting scheme including pionic effects. The scheme treats the 
pion mass and the inverse scattering length over the intrinsic short distance scale as small 
parameters. Working at next-to-leading order in this scheme, we show that the expression 
obtained is identical to one obtained using the recently introduced PDS approach which is 
based on dimensional regularization with a novel subtraction scheme. This strongly supports 
the conjecture that the schemes are equivalent provided one works to the same order in the 
power counting. 
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Weinberg [Q pointed out that there is a probiem with using the standard techniques of chirai perturbation 
theory (yPT) in nuciear physics (as opposed to hadronic physics). The difficuity is that nuciear physics 
has iow energy scaies quite distinct from m^^.. In iow energy s-wave scattering, the scattering iength, a, 
is extraordinariiy iong on hadronic scaies (~ 23. fm in the ^S'o channei and ~ 6 fm in the where 

m~^ ~ .25 fms). These unnaturaliy iong scaies set the scaie at which a pureiy perturbative treatment breaks 
down. This suggests that a conventionai perturbative approach breaks down at an absurdiy iow momentum, 
k ^ 1/a « 10 MeV in the ^S'o channel. Moreover, 1/a plays the role of the short distance scale (which we 
denote as A) in the perturbative problem so m,n. expansion will not converge. 

There has been considerable effort dealing with this “unnatural scattering length” problem during the 
past several years The goal is to implement the general underlying philosophy of effective field 

theory (EFT) to a problem in which m-,^ is not the only light scale. Weinberg’s proposal was to use 
yPT to develop a potential (two-particle irreducible graphs) which would be iterated to all orders yielding a 
Schrodinger equation. In order to implement chiral power counting for the two-particle irreducible graphs, 
four-nucleon contact interactions are required. However, as noted by a number of authors [pp|Jll|] when 
iterated to all orders, there is no straightforward way to implement the nonperturbative renormalization in 
a fashion consistent with the power counting. 

Two main approaches have emerged to get around this problem. The first is via cutoff EFT 
The scheme is to put in ad hoc short-ranged form factors in front of the contact interactions. By fitting the 
coefficients to get certain low energy observables (such as the scattering length) correctly one renormalizes the 
interactions in a way which is insensitive to the short distance physics. One can test this by demonstrating 
explicitly the insensitivity of observables to the cutoff scale [|^. The scheme, as implemented to date, has 
the disadvantage, however, that the calculation of observables has not been formulated directly in terms of 
an expansion parameter. This obscures one of the main advantages of EFTs, the ability to make a priori 
estimates of accuracy. 

The second approach is the so-called pole divergent subtraction (PDS) scheme of Kaplan, Savage and 
Wise j^. The scheme is based on an effective Lagrangian containing nucleons and pions and including 
nucleon contact interactions. In PDS the divergent integrals which arise from these interactions are treated 
via dimensional regularization with an unusual subtraction scheme: one subtracts poles at e = 1 from the 
analytically continued integrals rather than only poles at e = 0 as in the usual APS scheme. While this 
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prescription seems rather arbitrary, the consequence of this scheme is that coefficients in the Lagrangian 
have a systematic power counting which is not destroyed by large values of a, and that observables can be 
systematically expanded. A renormalization group analysis of the coefficients within this scheme shows that 
the approach requires that the lowest order s-wave interaction needs to be iterated to all orders—building 
in a nonperturbative treatment suitable for large scattering lengths—while all higher terms and all pionic 
effects are treated perturbatively. The principal advantage of this scheme is quite clear, one ends up with a 
treatment which is systematic in Q/A where Q k, rriT^ and thus in principle one can make a priori estimates 
of errors due to higher order effects. A number of observables have been calculated with this approach at 
next-to-leading order including nucleon-nucleon scattering [^ , and deuteron electromagnetic form factors 
pl| and polarizabilities |^2|. 

The purpose of this paper is to show that systematic power counting rules for observables can be imple¬ 
mented straightforwardly in the cutoff approach and to demonstrate the equivalence of the cutoff treatment 
with the PDS approach (provided both treatments implement the same power counting rules and work to 
the same order). The explicit demonstration of this equivalence is useful, if for no other reason, since it 
helps establish the validity of the formalisms. We explicitly demonstrate this equivalence in the context of 
scattering in the channel at next-to-leading order including the pionic effects. 

We implement a cutoff in configuration space by introducing a separation length scale, i?, which essentially 
serves as a renormalization scale. We parameterize the effects of the physics at lengths shorter than R and 
calculate the effects for lengths greater than R. The coefficients parameterizing the short distance physics 
vary with R in such a way as to offset changes in the long distance physics. In this way the formulation is 
renormalization group invariant. As a practical matter, R must be chosen to be sufficiently large that the 
effects of the short distance physics are contained almost completely within R. The scale of the long distance 
physics is and we impose a chiral expansion. Accordingly we must assume formally that rriT^R <C 1. 

We work with a chiral expansion with held fixed, which allows for the unnaturally large value of 
a. To simplify our power counting, we introduce a single low momentum scale, Q, and define the following 
power counting rules: 

m-n ^ Q Ija^Q k^Q (1) 


where k is the relative momentum of the colliding nucleons. This analysis based on power counting with a 
single low scale is consistent with the PDS analysis of ref. [Q. A couple of innocent assumptions go into the 
derivation of the expression for the phase shifts: i) at large distances the nucleon-nucleon interaction can be 
described by a Schrodinger equation ; ii) the energy dependence of the logarithmic derivative of the wave 
function at R is generically set by some short distance scale A. 

Weinberg’s power counting 0 at the level of the two-particle irreducible amplitudes is valid. We use the 
Fourier transforms of these amplitudes into configuration space. The implementation of the Q expansion is 
quite simple. One works in the Weinberg expansion up to some fixed order, including only the long distance 
pion exchange terms which we denote as Idong- By hypothesis R is taken to be sufficiently large that the 
effects of the short distance potential are essentially contained within R with an accuracy comparable to or 
better than the order to which we are working. Thus by assumption, for r > R, the entire potential is given 
by Mong- The associated Schrodinger equation is subject to boundary conditions at R, which come from 
fixing the energy dependence of the short distance physics. 

For simplicity, we consider the singlet channel and calculate kcot{6) at next-to-leading order—which is 
0{Q'^). The long range potential contributing at lowest nontrivial order is simply the central one-pion 
exchange potential which is given in configuration space by 


Idong(^) — h(,pep(^) 


-9a^I 

167r/2 


g-m^r 


r 


( 2 ) 


Note in this work we are using the convention that /^r ~ 93 MeV. Thus, the /tt used in this work is 1/^/2 
times /tt used in ref. |^. For r > R the system is described by the following radial Schrodinger equation 
for the relative position. 




-I- MViongir) - k^ u{r) = 0 


( 3 ) 
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where the s-wave component of the wave function is given by ips = w/r, and M is the mass of the nucleon. 

It is to useful express u{r) as an expansion u{r) = uo(r) + ui(r) + U 2 (r) + ... with uo(r) a solution of 
the noninteracting Schrodinger equation and 

Ui+i(r) = - MViong(r)ui(r) (4) 

It is immediately apparent that if the series for u converges, then any set of u^s which satisfy eq. (^) sums 
to a M which solves the original Schrodinger equation. The expansion is essentially a Born series except that 
it is only valid for r > i? and has nontrivial boundary conditions at r = i?. We impose Ui{R) = 0 for * > 1 
which puts all short distance physics effects into uq. Since uq is a solution of the free Schrodinger equation 
it can be written as uo(r') = sin(A:r + do); do simply represents what the singlet phase shift would have been 
had the long distance potential been removed. 

Simple dimensional analysis implies that each term in the expansion for u is order raT^jh. {QjK) smaller 
than the previous term; one expects that as —*■ 0, the expansion for u to converge. The leading nontrivial 
effect of the pions is contained in ui. It is straightforward to compute ui using Green’s function methods—the 
appropriate Green’s function is sin(fc(r — r'))/k so that 


u,(r) = M (5) 

We use the wave functions asymptotically in the determination of the phase shifts. Inserting uq = sin(/cr + do) 
into eq. (|) yields 

lim Mi(r) = sin(A:r) [cos(do)/i(A:) + sin(do)/ 2 (fc)] - cos(fcr) [sin(do)/i(A:) + cos(do)/ 3 (fc)] ( 6 ) 

r—^■oo 

where 

M 

fiik) = -^ J dr'Viong{r')sm{kr') cos{kr') 

M 

f2{k) = -^ J dr'Viong{.r')cos^(kr') 

M r°° 

fsik) = -^ J dr'Viong{r')sm'^{kr') 

Asymptotically the wave function is of the form Asin(A:r + d) where A is the overall normalization. From 
this asymptotic form and the fact that to this order u = uq + ui one can obtain an expression for fccot(d): 


( 7 ) 

( 8 ) 
(9) 


kcot{6) = 


fccot(do) (1 + /i) + fc /2 
1 - /i - /ccot(do) ifa/k) 


( 10 ) 


We need to establish how the phase shift relation above behaves for small Q. First, we note that we can 
expand /ccot(do) as a function of k^. This is an effective range expansion for a system consisting only of the 
short distance part of the interaction: 

fccot(do) = -l/oshort + 1/2 + ^2 + U 3 /c® + /c® + ... ( 11 ) 

By hypothesis, the scale of the coefficients are generically set by the short distance scale A: r® ~ A“^, 
u/ ~ This reflects the assumption that at the matching distance, R, the energy dependence of the 

logarithmic derivative of the wave function is set by A. Generically, one would expect l/oghort ^ A but here 
we will allow for the possibility that l/oshort is unnaturally small and will assign it a power counting scale 
of Q. As we will see subsequently 1/a = l/oshort + 0{Q'^/A). We are working to order in /ccot(d) so we 
truncate the expansion at the effective range term. It is useful to rewrite 1 /aghort in terms of its value in the 
chiral limit plus a chiral correction 
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( 12 ) 


I 7 2 

- =-h dm^ 

^short ^0 

where we have truncated the chiral expansion at order 

Nontrivial chiral behavior is contained in the / functions. Up to an overall multiplicative factor of 
{g\M)/{IQ tt f"^) the three /’s are functions k with parameters m,r and R and mass dimension zero. One can 
write such functions as 


f^{k;m^,R) ~ m, 


-AM 

167r/2 


f^{k/mT,-m^R) 


(13) 


The preceding form is convenient in that we can do Q power counting treating 771,7 as a small parameter 
of order Q/A. Accordingly we expand in to,,!?; the lowest order contributions come from taking m-^R — 0. 
Given a fixed to,,, this amounts to taking R to zero in the integrals. For /i and /s this can be done quite 
straightforwardly: 

/i(fc) = -^ J dr'Uiong(r')sin(fcr') cos(fc/) (l + e>(TO„i?)) 


g\M fm 




M /■°° 

fsik) = — J dr'Viong{r')sm^{kr'){l + 0{m„R)) 


= -rrij. 


g\M /TO„\ , 


64 -Kfl 


(x) 


In 1 + 


(l + C>(TO„i?)) 


(14) 


The corrections due to finite TO„i? contribute to the final results only at higher order in Q and will be 
dropped here. In evaluating /2 one finds a logarithmic divergence when taking R to zero. It is useful to 
rewrite the cos^ in /2 as 1 — sin^ so as to re-express /2 as 


M r°° 

f 2 {k) = -h{k) + -r dr'Uiong(r') 

R 


= -fsik) + rriT, (^) ln(TO„i?)) (1 -b 0{m^R)) 


16 tt/, 


(15) 


The logarithmic dependence on R in this expression will ultimately be absorbed in the renormalization 
process—by renormalization group invariance d has a ln(i?) dependence against which the dependence in /2 
cancels. 

We can now insert our expressions for A:cot(i5o) and the fi functions into eq. ( p^ and expand the resulting 
expression in Q up to to obtain: 


A:cot(<5) = —— + 
ao 


ml (d -b 7 + ln(TO„i?)) + - i2 


1 glM 
A MTTf^ 


k^ 


In ( 1 H-^ 

mi 


qlM /TO„\ , f 2k\ , 

I ^ \ / .] + mi In 1 + — 


— tan“^ (- 

On 167r t£ \ k y \ TO,, , 


AM 


mi 


ao 167r/2 \ k J \m-^ J M-nfl 

Note that from eq. & and the fact that fccot(d) |^,^q = —1/a one finds that 


1 

a 


1 

ao 


2 9aM 
" 1677/2 


(d -b 7 -b ln(TO„i?)) -b 


g\M f 2to„ 1 


16^/2 


ao 


1 

ao 


- - = - + o(Q 7 A) 


(16) 


(17) 


We wish to compare our result with that of the PDS scheme of Kaplan, Savage and Wise Q. We will 
not review the scheme in detail here since ref. [Q does a good job of explaining clearly the steps necessary 
to implement the scheme. 
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We take the relevant expressions from ref. ||^ and use the notation of ref. (except we convert to our 
normalization of /^). In the PDS scheme, one calculates the scattering amplitudes directly order by order 
in Q, and infers the phase shifts or kcoi{5) through them. By unitarity, the general form of the amplitude 
is given by 


47r 1 

M (k cot(i5) — ip) 


(18) 


In the ^So channel, apart from pion exchange there are three nucleon contact operators in the Lagrangian 
which contribute at this order—there is contact interaction with no derivatives, a chiral correction to it and 
two derivative interaction, parameterized by (scale-dependent) coefhcients Cq, D 2 and C 2 , respectively. The 
amplitude at order Q~^ is obtained by iterating the Cq contact term to all orders. 


■4-1 — 


-Cn 


47r 


1 


[l + Q^{ik + p)] Mil/ao + ik) 


(19) 


where the second form can be taken as a definition of uq. This definition is consistent with the definition 
used in the cutoff approach; both correspond to the value of the nucleon-nucleon scattering length in an 
artificial world where the quark masses were strictly zero. It is also useful at this point to observe that the 
PDS formalism requires that I/uq is formally treated as being of order Q since A-i is order 1/Q. At order 
the scattering amplitude is given as the sum of five terms: Aq = Aq'^ + -I- -I- A^^^ A^^\ 

They are given by 


A\!'^ = -C 2 fc" 


- 

v^n — 


[A. 


1 2 


Cn 


A 


(III) _ 


a! 


{IV) 


a: 


iv) 



ik) 

TTItt 2fc 


p + ik 


tan 


i tan 



= -Dom: 


where /i is the renormalization scale. Equating 1 over the scattering amplitude from eq. (^ 
Ao) and expanding to order gives |M 


( 20 ) 

with l/[A-i + 


kcot{6) = ik 


47r 


47rAo 


MA-i MAi^ 

Inserting the expressions for A-i and Ao into eq. (pi]) and simplifying the resulting expression yields: 


( 21 ) 


A:cot(i5) =-h 

Oo 


9aM 


mi In 




- + 


1 
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9\M 
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1 g\M 
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1 -k 
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4fc2 
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( 22 ) 


where to aid comparison with the eq. di) we have used eq. dil) to eliminate Co in favor of ao- Comparing 
eq. (0) with eq. (p^ we see the two are identical provided that we identify the quantities in square and 
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curly brackets in eq. with (g‘\M)/(IGirf^) (d + 7 + ln(TO^i?)) and rg/2, respectively. The quantities 
in square and curly brackets in eq. ( P2[) contain coefficients from the Lagrangian which may be adjusted to 
reproduce the correct short distance physics and hence may be matched to the analogous terms in eq. (|^). 
This completes the demonstration that expressions for ^ So scattering derived using the Schrodinger equation 
in configuration space is equivalent to the PDS result. 

The analysis based on the configuration space Schrodinger equation is along the line of conventional poten¬ 
tial model treatments. We simply imposed a power counting scheme in Q on the traditional approach. The 
PDS approach is superficially very different. Our demonstration of the equivalence confirms the supposition 
that the PDS approach describes the same physics provided one is in a regime in which the power counting 
in Q is valid. The extent to which we are in such a regime will be addressed in a subsequent paper. 
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acknowledges the support of the U.S. Department of Energy under grant no. DE-FG02-93ER-40762. 
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